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In modeling low-dimensional electronic nanostructures, the evaluation of the electron-electron 
interaction is a challenging task. Here we present an accurate and practical density- functional 
approach to the two-dimensional many-electron problem. In particular, we show that spin-density 
functionals in the class of meta-generalized-gradient approximations can be greatly simplified by 
reducing the explicit dependence on the Kohn-Sham orbitals to the dependence on the electron spin 
density and its spatial derivatives. Tests on various quantum-dot systems show that the overall 
accuracy is well preserved, if not even improved, by the modifications. 
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I. INTRODUCTION 

With the present technology the electron gas can be 
confined in various ways to create nanoscale devices of 
lower dimension. The field of two-dimensional (2D) 
physics has grown rapidly alongside the development of 
electronic devices such as quantum Hall bars and point 
contacts, and semiconductor quantum dots. When mod- 
cling these systems, the finite extent in the growth direc- 
tion (say z) can often be neglected, so that the system 
is well described by a 2D Hamiltonian in the effective- 
mass approximation^ In this respect, the building block 
is the 2D electron gas (2DEG), whose properties are well 
known in the literature^ 

Density-functional theory (DFT) and its extensions 
have become the method of choice to describe the 
electronic properties of three-dimensional (3D) systems 
such as atoms, molecules and solids 4^ Despite the fact 
that the many 3D density functionals developed for the 
exchange-correlation (xc) energy and potential fail in 
the quasi-2D limit^— the derivation of explicitly 2D xc 
functionals has started only very recently^— In finite 
2D systems, most of these functionals overperform the 
2D local spin-density approximation (LSDA), which is 
a combination of the analytic exchange energy of the 
2DEG— and the corresponding correlation energy hav- 
ing a parametrized fori m 19 : 20 The encouraging results 
obtained with the new functionals indicate that DFT in 
2D has entered in a more mature phase. 

Among the newly proposed functionals, our focus 
is on those approximations that have as ingredients 
the electron density and its spatial derivatives, the ki- 
netic energy density and the paramagnetic current^— 
In other words, the expressions of these functionals 
are current-dependent meta-generalized-gradient approx- 
imations (meta-GGAs), and therefore they explicitly de- 
pend on the Kohn-Sham (KS) orbitals. As it is valid 
for 3D systems, also for 2D systems the meta-GGAs are 
very accurate. However, the price for their accuracy is 
the more involved numerical implementation (if applied 



self-consistently) as well as the (case-dependent) numer- 
ical burden in the applications. In order to simplify both 
the mentioned tasks, we explore here to which extent and 
how the explicit dependence on the KS orbitals may be 
reduced to the dependence on the electron density and 
its spatial derivatives, yet possibly maintaining a satis- 
factory level of accuracy. In other words, we examine the 
path from implicit to explicit density functionals in the 
class of meta-GGAs. A similar study has been already 
carried out for 3D systems, introducing some Laplacian- 
level meta-GGA functionals. 21 With the present work, 
we explore this possibility for low-dimensional systems. 
We find that the performance after the simplifications is 
well preserved, and in some cases even improved. 



II. REVIEW AND MODIFICATIONS OF THE 
FUNCTIONALS 



In the following, we review the main ingredients of 
recently derived functionals and suggest how their ex- 
pressions may be significantly simplified in a consistent 
fashion. Ideally, the ultimate goal is to obtain the best 
performance with the least (numerical) effort. In par- 
ticular, we consider (i) an exchange-energy functional 
obtained through the modeling of the exchange-hole (x- 
hole) functions? 9 - (ii) an exchange-energy functional ob- 
tained from the one-body-density-matrbsii (1BSDM), 
and (iii) correlation-energy functionals obtained through 
the modeling of the correlation hole i 12 ' 13 Details of the 
derivation of the functionals can be found the in the men- 
tioned literature. 
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Exchange energy from the exchange-hole 
functions 



The exchange energy E x can be expressed through the 
x-hole functions as^ h%(rx,T2) 

•y p />00 />27T 

E -=2^J * r P°W J Q ds J <^>K(r,T+s), (1) 



with 



(2) 



where V'fco-M are the KS (spin) orbitals. From Eq. JT]), it 
is apparent that an approximation of the x-hole also pro- 
vides an approximation for the exchange energies. Equa- 
tion ([T]) also suggests that the details of the angular de- 
pendence of the x-hole are energetically negligible: all we 
need is the cylindrical average of the x-hole, i.e. h x (r; s). 
As the basis of our model^ we have chosen 
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K(r;s) = — I d<j> s Z£(r;r + s) 



--exp [-a(r) (b(r) + s 2 )] 
Io (2a(r)VKrTs) , 



(3) 



where Iq(x) is the zeroth-order modified Bessel func- 
tion of the first kind. This model provides the cor- 
rect sign of the x-hole and the correct normalization 
J d 2 s h%(v, s) = — 1. The non-negative functions a(r) 
and 6(r) are introduced to reproduce the short-range be- 
havior of the x-hole. It is important to note that the 
curvature of the x-hole in 2D is given by 9 - 



CJ(r) = \ 



2 Ax(r) 



P<x( r ) 

where r CT is (twice) the spin-dependent kinetic-energy 
density, and j Pi(T is the spin-dependent paramagnetic cur- 
rent density. 

In Ref. |3, we have applied the above scheme in two 
ways. In the first instance, we have employed it in its full 
spirit by numerically determining a and b at each point 
in space. Then, in a fully self-consistent application, 
the scheme should be implemented within the optimized- 
effective-potential (OEP) method.2 2 - 2 ^ This is necessary 
because r CT and j PiCr (in non-current-spin-density func- 
tional calculation s 27 ' 28 ) make the overall scheme explic- 
itly dependent on the KS orbitals. As a consequence, the 
corresponding numerical task would be highly nontrivial. 

In the second instance, we have analyzed the 2DEG 
limit, for which 



(5) 



In this way, a local density functional has been recovered, 
which was seen to improve over the exchange energies 



obtained within the standard LSDA when applied to few- 
electron quantum dots. 

Next we proceed from the review of the functional to its 
modification. As discussed above, the task is to remove 
the explicit reference to the KS orbitals, yet retaining 
some degree of flexibility in dealing with inhomogeneous 
systems. This may be achieved by introducing the fol- 
lowing modification ! 15 ' 29 ! 30 

N " 1 i 2 (r) 

r CT (r) =5>fc ff (r) v(r) = 2^p 2 (r) + -V 2 p CT (r)+ J -^ 

k=i 6 pa[r) 

(6) 

Correspondingly, we obtain 



C£(r) C£(r) = -Tr^(r) + ^V 2 P(T (r) 



1 (Vp.(r)) 1 



8 Pa(r) 



(7) 



It is worth mentioning that due to the last (current 
dependent) term on the right hand side of Eq. ([6]), the 
modified x-hole curvature manifestly preserves its gauge 
invariance. As described in Ref. [^, the functions a and b 
are determined from 



a = npa- exp(y), 



and 



y 



ex p (-y) 



(8) 



(9) 



where y = ab satisfies 
CZ 



(y - l)exp(y) 



-1 



1 V 2 Pa , 1 (Vp) 2 



12tt 



P'n 



Sir p\ 

(10) 

If a solution does not exist, we set£ y = 0. This corre- 
sponds to the 2DEG limit mentioned just above. From 
Eq. (|T0|) . it is apparent that the second and third terms 
are relevant when the curvature and the gradient of the 
spin-density are non- negligible. 

Going back to the first modification r CT — > f a , we 
refer to the new (angular-averaged) x-hole function as 

h%{r,s). The corresponding x-hole potentials denoted as 

^x.model( r )> read as follow^ 



U° 



odclV 



(r) = 27r/ dsK(r,s) 



from which, the exchange energy is obtained as 



El 



odcl 



d 2 rp a {v)U:(r) 



(11) 



(12) 



The last two quantities calculated without the modifica- 
tion of r a are denoted below simply without the tilde 
symbols, i.e., £^ modcl and £% model . 

Finally, we point out that the above modifications of- 
fer a straightforward way to calculate the correspond- 
ing KS exchange potential functional derivative, 
u£ = 5E x /{8p a ). The properties and performance of 
these potentials will be assessed elsewhere. 
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B. Exchange energy from the 
one-body-spin-density matrix 

Another way to express the exchange energy is to make 
use of the 1BSDM j a (ri,r 2 ), 3 ' 4 so that 



E ,. 



-E 

2 ^ 

2tt 



dfr 



ds 



s 

2' r 



with 



7<r(ri,r 2 ) 



E 

fc=i 



(13) 



(14) 



where V'fecr( r ) are the KS (spin) orbitals. Clearly, an ap- 
proximation for the 1BSDM implies an approximation 
for the exchange energy. Also, it is apparent that the an- 
gular dependence of the 1BSDM is energetically of minor 
importance. Therefore, as a basis of our approximation 
we have considered the following expression^ 



1 

2^ 



&(r) 



f S 

v r+ 2' r -2, 



A a {N a ) 



(15) 



where j3(r) is chosen to reproduce the exact-short behav- 
ior of the 1BSDM 



lr CT (r) lV^(r) 1/j 
2 



2p CT (r) 8 p a (r) 



/v(r) 



(16) 



and A a (N a ) is obtained through the normalization of the 
particle number for each spin channel. In Rcf. El we have 
employed this scheme leading to accurate results for vari- 
ous quantum-dot systems. In addition, we have observed 
that this approach and the one of Sec. Ill Al coincide in 
the 2DEG limit. 

Here we suggest to simplify the present functional by 
making use of Eq. © in Eq. ([TBI . This yields 



p- 1 (r) = np a (r) 



1 W(r) 
24 PrT (r) 



(17) 



We emphasize that no gradients of the spin-density ap- 
pear in this expression. The resulting expression is man- 
ifestly gauge invariant. The relevance of the gauge in- 
variance of the expression in Eq. (|16p has been already 
verified in Ref. !ll). The final expression reads as 



E 



1BSDM 



d 2 - 



= -\ E ' ,r " 

£T=t4 



(18) 



where A a is determined through the normalization as 



N a = n 



d 2 r 



l + 2A 



(19) 



Equation (1181) together with Eq. (1191) provide another 
density functional for the exchange energy. Finally, the 
x-hole potential has a form 



z,lBSDM 



(r) 



r- 3 I 
+ 4 ° 



P Ar)fc /2 (r). (20) 



Similarly to the previous section, the exchange-hole po- 
tentials and exchange energies calculated without the 
modification of t ct are denoted below simply without the 
tilde symbols, i.e., l/' 1BgDM and ££ 1BSDM . 



C. Correlation energy from the correlation-hole 
functions 



High predictive power in the application of DFT re- 
quires the accurate treatment of the electronic correlation 
in both inhomogeneous systems and in the limit of the 
homogeneous electron gas. We have achieved this goal 
in 2D by generalizing our previous approximation^ to a 
parameter-free form^ which reproduces the correlation 
energy of the 2DEG while preserving the ability to deal 
with inhomogeneous systems (quantum dots). 

The correlation energy is expressed in terms of the 
cylindrical average of the (coupling-constant dependent) 
correlation-hole (c-hole) functions h a c \ (r, s) as follows: 3 ' 4 

//>oo />1 
dv Pry (v) / ds / d\hg(T,s) . (21) 
Jo Jo 

It is obvious that an approximation for /i™(r, s) implies 

an approximation for E^ a . In modeling these quantities, 
we have proposed a formic 



(s - z aa (v)) D a {r) 



1 + | A2w(r) 
s - z CTff (r) 



exp ( 


9irs 2 \ 







(22) 



1 + 2\z< r9 {r) 



exp 



44 ff (r) 



(23) 

for the same- and opposite-spin cases, o~cr' = era and 
ao-' = aa, respectively. Here 



2 \ 4 p a p a 



(24) 



and 



z aa {v) := 2<w|[/*(r)|-\ (25) 
z aS {v) := c^[\K{v)\- l + \Ul{v)\- 1 ] . (26) 

Equations (|2"5j) and (12"B"|) are proportionality relations that 
may be enforced locally in space [see Eqs. (j27[) and 
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below]. It is apparent that z aa i{y) set the characteristic 
sizes of the c-hole functions in terms of the sizes of the 
x-hole functions. The idea behind this assumption is the 
following: the smaller the x-hole around each electron is, 
the more tightly the electrons are screened. Therefore, 
they are expected to be correlated much less. 

The above modeling provides 12 : (i) zero correlation en- 
ergy for one-particle systems (as the exact one); (ii) exact 
short-range behavior of the A-dependent c-hole functions; 
(iii) a "reasonable" decay in the limit s — > oo; (iv) exact 
normalization of the A-dependent c-hole functions. Fur- 
thermore, (v) c aa i can be defined in such a way that the 
total correlation energy of the 2DEG is exactly repro- 
duced.— As a result, when the (average) density has a 
realistic range, < r s = < 20, we can use the 

following (approximate) parameterizations: 



a log(r s ) + fir 



7 



(27) 



with a = -0.1415 1, /3 = 1.226 1, 7 = 0.144 99, and 

caa[r s ] = Sri (28) 

with 5 = 0.663 25 and £ = 0.123 96. When using the 
present correlation functional, the coefficients c^o-'^sKr) 
must be calculated at each point in space by making use 
of the density, that is, r s (r) = l/\Ar/9(r). 

Here we propose a few simplifications along the lines of 
the previous sections. First, we may replace U° with an 
approximate expression obtained in Sec. Ill Al Secondly, 
we apply Eq. (|6|) to Eq. (l24l) leading to 

D 9 {r) -> D.{t) = npl(r) + Vp ff (r) - J^^TT^ • 

6 8 pair) 

Now, conditions (i) and (ii) given above are no longer sat- 
isfied. The latter modification clearly affects the same- 
spin c-hole functions. Applying both of the described 
simplifications - which is naturally required in order to 
make the functional orbital-free - the correlation energy 
can be expressed in terms of the (spin-dependent) c-hole 
potentials, U° a (r), as follows: 



Er = 



1 



with 



16 



drp a (r)Ur (r) , 



370^(1-)^ (r) 



(30) 



817T 

2z aa (v) - 3 In ( |z CTCT (r) + 1 



(31) 



and 



WW = (2-<Mr) 

x [22^(r)-m(2iw(r) + l)] , (32) 

where z aa i (r) are obtained by replacing U° with U£ in 
Eqs. (H5J) and (EU). 



III. TESTING THE MODIFICATIONS 

Next we test the modifications for different 2D 
quantum-dot systems. As a standard test set we con- 
sider parabolic (harmonic) dots consisting of N electrons 
confined in an external potential v ex t(r) = uj 2 r 2 /2. First 
we use the octopus code^i to solve the KS problem 
self-consistently using the exact-exchange (EXX) func- 
tional within the Krieger-Li-Iafrate (KLI) approxima- 
tion.— Then the resulting KS orbitals - and for the 
modified functionals solely the electron density and its 
gradients - are used to compute the energy expressions 
and their ingredients introduced in the previous sections. 
The obtained exchange energies can be directly compared 
with the EXX- KLI results E® xx . In addition to this test 
set, we also consider a large quantum dot where we com- 
pare with the LSDA, as well as a rectangular quantum 
slab. In the case of correlation, we exploit the numeri- 
cally exact configuration-interaction data^ for the total 
energies ££ot > so that the reference correlation energy can 



be evaluated from El = E, 



:ef 



E, 



EXX 



A. Exchange energies 

First we test the ingredients entering in the expressions 
for the exchange energy in Sees. Ill Al and ITTBl 

Figure [TJa) shows the original and modified kinetic- 
energy densities [as defined in Eq. ©] of a spin-polarized 
three-electron parabolic quantum dot with uj = 1/4. The 
characteristic step in r CT at the shell of the quantum dot 
at r ~ 3 is significantly smoother in f CT . However, this 
difference is partly washed away in the local curvature of 
the x-hole shown in Fig. [ljb). 

In the exchange-hole potential shown in Fig. (Ha) - 
computed with the functional described in Sec. IIIBI - 
the difference is reduced further, so that the results are 
almost identical. Moreover, they agree very well with 
the EXX-KLI result corresponding to the Slater potential 
(dotted line). Naturally, this similarity leads to precise 
exchange energies as explicitly shown below. 

In Fig. [5{b) we show the exchange-hole potentials cal- 
culated with the functional described in Sec. Ill Bl Again, 
the results with and without the modification of the 
functional are similar, although the relative difference is 
larger than in the previous case [Fig. HJa)]. Both poten- 
tials, however, deviate rather strongly from the EXX-KLI 
result. This tendency is already present in the original 
functional, which has been tailored mainly to produce 
accurate exchange energies^ which is indeed the case as 
shown below. 

Table U shows the exchange energies calculated for sev- 
eral quantum dots. The set includes four cases with or- 
bital currents (rows 1, 3, 4, and 7), in two cases aris- 
ing from an external magnetic field perpendicular to the 
2D plane (rows 3 and 7)£& We consider modifications 
for both £™ odcl and £^ BSDM , respectively. Overall, the 
modifications preserve the excellent performance of the 
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FIG. 1: (color online) (a) Comparison of the original (solid 
line) and modified (dashed line) kinetic-energy density of a 
spin-polarized three-electron parabolic (cj = 1/4) quantum 
dot. (b) Local curvature of the exchange hole calculated from 
the original (solid line) and modified (dashed line) kinetic- 
energy density. 



functionals very well (see the last row of Tabled . For the 
1BSDM approximation the modification even improves 
the performance. The LSDA is giving clearly the worst 
accuracy of the tested functionals. 

In addition to the test set of Table Q] that covers only 
few-electron quantum dots, we now consider two rather 
different cases. First we focus on a large 48-electron 
parabolic quantum dot with uj — 0.3373 at a magnetic 
field of B = 3.05 T. This partially spin-polarized (total 
spin S = 3) ground state has a compact "spin droplet" 
on the second-lowest Landau level, and its existence has 
been confirmed in recent spin-blockade experiments 
Here we have performed a LSDA calculation and use 
that density as an input in the functionals. Figure [3] 
shows the kinetic-energy densities and exchange-hole po- 
tentials calculated with the functionals obtained from the 
modeling of the exchange hole (see Sec. Ill A|) . The orig- 
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FIG. 2: (color online) Exchange-hole potentials of a spin- 
polarized three-electron parabolic (uj — l/4)quantum dot. (a) 
Result of the functionals in Ref.0 without (solid line) and with 
(dashed line) the modification (se e Sec. Ill A")) , (b) The same 
as (a) for the functional in Ref . ITU (see Sec. Ill B] ). The dotted 
line corresponds to the EXX-KLI result. 

inal and modified T a for both spin-up and spin-down 
electrons are very similar, and the resulting exchange- 
hole potentials are practically the same. The exchange 
energies are E™ odcl = -22.18 and £™ odcl = -22.25 
(difference of 0.3%). In comparison, the LSDA yields 
^LSDA = —21.11. In lack of a reliable EXX reference 
result for a system of this size it is not possible to judge 
whether the exchange energy from the model(s) or from 
the LSDA is more accurate. However, knowing that the 
LSDA typically underestimates the (absolute value of) 
E x , our results for Ef odcX and £™odei deviate from the 
LSDA in the correct direction. Most importantly, the 
exchange-hole potentials from the model are more accu- 
rate, especially in the asymptotic region. 9 

Our second example of a quantum-dot system that 
differs from those in Table U is a 16-electron rectan- 
gular hard- wall quantum slab with size 2\Z2tt x v2~7T 
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TABLE I: Exchange energies for fully spin-polarized parabolic 
quantum dots calculated using the functional of Ref. 
{E™ odcl ), its modification described in Sec. [2] (E7° dcl ), the 
functional in of Ref. El (£ 1BSDM ), its modification (i^ BSDM ) 
described in Sec. Ill Bl and the local spin-density approxima- 
tion (^ SDA . They are compared with the EXX-KLI result 
(_E EXX ), so that the last row shows the mean absolute error 
in percentage. 



N 


B(T) 


j-rmodcl 
&x 


jjimodcl 
&x 


^jlBSDM 


&x 


j-iLSDA 
E x 


pEXX 
&x 


2 





-0.626 


-0.634 


-0.618 


-0.620 


-0.583 


-0.626 


3 





-1.038 


-1.043 


-1.029 


-1.021 


-0.963 


-1.021 


3 


2 


-1.038 


-1.056 


-1.029 


-1.037 


-0.979 


-1.039 


4 





-1.421 


-1.435 


-1.416 


-1.408 


-1.332 


-1.374 


5 





-1.865 


-1.876 


-1.846 


-1.842 


-1.745 


-1.816 


6 





-2.267 


-2.275 


-2.249 


-2.241 


-2.126 


-2.214 


6 


3 


-2.349 


-2.391 


-2.344 


-2.362 


-2.241 


-2.357 






1.5 


2.4 


1.4 


0.9 


4.9 





corresponding to ~ 90 nm x 45 nm in SI units.— In 
Fig. 2] we compare the exchange- hole potentials given by 
the 1BSDM functional (see Sec. HTB|) to the EXX-KLI 
(Slater) potential computed. The overall shapes are very 
similar, but as expected, the EXX-KLI potential is con- 
siderably smoother. Interestingly, however, the modified 
potential is qualitatively closer to the EXX-KLI result 
than the original one. Regarding the exchange energies 
both functionals perform similarly: iji BSDM = 13.13, 
^ibsdm = 13 15) and e exx = 12 7 ThuSj the modi _ 

fixation in r a is well justified also when considering a 2D 
system with a hard- wall geometry. 



B. Correlation energies 

Figure [5] shows the spin-pair components of the 
correlation-hole potentials for a spin-unpolarized (total 
spin S = 0) six-electron parabolic quantum dot (lu — 
1/4) calculated with the model of Ref. [TH in comparison 
with its modifications introduced in Sec. Ill CI The mod- 
ification induces clear devitations from the original po- 
tential, especially for the same-spin component affected 
by modifications in D a . For example, the bump at r ~ 4 
is due to the change of sign in D a in that regime. In con- 
trast, the opposite-spin component is independent of D a 
[see Eq. (f3"2"|) ]. so that the modified functionals are almost 
the same; here the choice of f^ mo( jei instead of C^exx 
has a negligible effect (solid and dashed lines overlap). 
In lack of an exact reference result we cannot assess the 
quality of the correlation- hole potential(s). Hence, in 
the following we will focus on the correlation energies 
for which reference results can be obtained as described 
in the beginning of Sec. IIIII 

In Table |TT] we test the effect of the modifications for 
the accuracy of the correlation-energy functional. As dis- 
cussed at the end of Sec. Ill CI we consider two approxi- 
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FIG. 3: (color online) (a) Original (solid lines) and modi- 
fied (dashed lines) kinetic-energy densities for spin-up and 
spin-down electrons in a 48-electron parabolic (u — 0.3373) 
quantum dot at B = 3.05 T. The dotted lines show the spin 
densities, (b) Resulting spin-up and spin-down exchange-hole 
potentials using the functionals obtained from the modeling 
of the exchange hole (see Sec. Ill A| . 



mations, where we either do not approximate U x in the 
same framework but use the EXX result, or then we apply 
the modification also to U x . Interestingly, the best result 
- apart from the original functional which is very close 
in accuracy - is given by E C (U£) for both spin-polarized 
and unpolarized cases. This finding may demonstrate 
the compatibility between the corresponding exchange- 
and correlation-energy functionals. Nevertheless, all the 
functionals introduced here are superior to the LSDA, 
whose error is an order of magnitude larger] 12 ' 13 

Finally, in Fig. [5] we plot the relative errors of the 
correlation-energy functionals as a function of N (cf. Fig. 
2 in Ref. fl3l) . It is interesting to note that, at least for this 
set of systems, the modified functionals show consistent 
behavior as a function of the number of electrons. Thus, 
it may be expected that the good performance continues 
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FIG. 4: (color online) Exchange-hole potentials of a 16- 
electron rectangular quantum dot calculated with the original 
(a) and modified (b) functionals obtained from the one-body 
spin-density matrix (see Sec. Ill B|l in comparision with the 
exact-exchange (Slater) potential (c). 



further to larger N. Unfortunately, a throughout test- 
ing of this is beyond the capability of numerically exact 
methods to provide accurate reference data. 



IV. CONCLUSIONS AND OUTLOOK 

In this work we have explored the possibility to 
modify meta-generalized-gradient approximations (meta- 
GGAs) for the exchange and correlation energies of two- 
dimensional systems to Laplacian-level meta-GGA ones. 
We have analyzed the effects of the according modifica- 
tions on various systems. Although the differences in the 
kinetic-energy densities can be considerable, the func- 
tionals considered in this work preserve well the quality 
of the exchange- and correlation-hole potentials, and in 
particular the corresponding energies. Overall, we find 
that the performance is well preserved, if not even im- 
proved, by the modifications. 

The simplified meta-GGAs provide significant prac- 
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FIG. 5: (color online) Spin-pair components of the 
correlation-hole potential for a spin-unpolarized six-electron 
parabolic (uj = 1/4) quantum dot calculated using the func- 
tional obtained from the correlation-hole modeling^ (solid 
line), its modification described in Sec. Ill Cl with LT^exx (ob- 
tained within the KLI approximation) as an input (dashed 
line), and with fX£ model as an input (dotted line). 



TABLE II: Correlation energies for spin-polarized (S = N/2) 
and unpolarized (S = 0) parabolic quantum dots calculated 
using the functional in Ref. [l3| (see Sec. Ill C|l , its modifica- 
tion with LT^exx (in the KLI approximation) as an input, the 
modified form with U£ as an input, and the local-spin-density 
approximation.— The last column shows the numerically ex- 
act reference result. The last row shows the mean absolute 
error in percentage. 



N 


S 


u 


E c 




-^c(f/x,modcl) 


jjLSDA 


E%* 


2 


1 


1/4 


-0.0115 


-0.0073 


-0.0085 


-0.0345 


-0.0100 


3 


3/2 


1/4 


-0.0225 


-0.0189 


-0.0200 


-0.0564 


-0.0226 


4 


2 


1/4 


-0.0399 


-0.0337 


-0.0330 


-0.0730 


-0.0337 


5 


5/2 


1/4 


-0.0570 


-0.0465 


-0.0468 


-0.0929 


-0.0484 


6 


3 


1/4 


-0.0681 


-0.0617 


-0.0613 


-0.1125 


-0.0640 


6 





1/4 


-0.390 


-0.379 


0.380 


-0.458 


-0.396 


3 


3/2 


1/16 


-0.0138 


-0.0106 


-0.0122 


-0.0382 


-0.0167 


5 


5/2 


1/16 


-0.0348 


-0.0278 


-0.0299 


-0.0659 


-0.0357 


6 


3 


1/16 


-0.0426 


-0.0372 


-0.0393 


-0.0796 


-0.0459 


6 





1/16 


-0.228 


-0.214 


-0.220 


-0.282 


-0.250 








9.5 


15 


11 


99 





ticality and numerical efficiency in the application of 
the functionals. Therefore, a self-consistent and multi- 
purpose implementation of the present toolbox of func- 
tionals is now within reach, enabling the investigation 
of (quasi-) two dimensional electronic nanostructures of 
experimental and technological relevance. 
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